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We prove that the amplitudes for the (d,p), (d,pn) and (e, e'p) reactions determining the asymp- 
totic behavior of the exact scattering wave functions in the corresponding channels are invariant 
under unitary correlation operators while the spectroscopic factors are not. Moreover, the exact 
reaction amplitudes are not parametrized in terms of the spectroscopic factors and cannot provide 
a tool to determine the spectroscopic factors. 
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Introduction. — Since the dawn of nuclear physics (d,p) 
reactions have been the main tool for extracting spec- 
troscopic factors (SFs) (we call them phenomenologi- 
cal SFs), which were compared with predictions of the 
independent-particle shell model (ISPM) (see [l| and ref- 
erences therein) . Later on electron-induced breakup reac- 
tions (see [2j and references therein) and nucleon knock- 
out reactions (see Q and references therein) became new 
tools to determine the SFs. Reduction of the phenomeno- 
logical SFs deduced from the (e, e' p) reactions compared 
to the IPSM ones has been discussed in 12J. Recently a 
similar reduction has been observed for the SFs deter- 
mined from the analysis of the single-nucleon knockout 
nuclear reactions Q as well. There are different possible 
sources of this reduction: the single-particle approxima- 
tion to the overlap function, antisymmetrization effects, 
an ambiguity of the optical potentials, accuracy of the 
DWBA, contribution of the coupled channels or the ef- 
fect of the short-range NN correlations. 

In the IPSM the model wave functions are given by the 
linear combination of the Slater determinant wave func- 
tions. However, the such model wave functions don't take 
into account short-range NN correlations caused by the 
repulsive core in the TV TV potential. In [J] the so-called 
unitary correlation operator method (UCOM) has been 
developed, which allows one to correct the IPSM taking 
into account the short-range repulsive NN correlations 
by applying unitary correlation operators (UCOs) onto 
the trial functions. In this work we consider the effect of 
such short-range correlations and demonstrate that the 
exact reaction amplitude is invariant under finite-range 
unitary transformation of the wave functions. 

(d, p) reaction amplitude. — We start our consideration 
from the deuteron stripping reaction A(d,p)B. Its am- 
plitude in the post-form is given by 



M dp =< X { f\B\V pB - U pB \^ ( + ] > 



(1) 



Here \1/,- ' is the exact d + A scattering wave function 
which has outgoing waves in all the open channels and 
satisfies 



Xf is the distorted wave in the exit p + B channel. It 

is a solution to {T pB + T B + U pB + V B )x ( f ) = E X ( } '\ 
where T pB is the kinetic energy operator for the relative 
motion of the proton and the center of mass of nucleus 
B while T B is the kinetic energy operator for internal 
motion of the nucleons in nucleus B. In addition, ip B 
is the nucleus B bound state wave function and Vy is 
the interaction potential of nuclei i and j, Uij is their 
optical potential, V B is the interaction potential of the 
nucleons in nucleus B. Note that the potential V pB may 
consist of the sum of the two- and three- nucleon forces. 
In the paper we will take into account only two-nucleon 



forces, i.e. V pB 
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the proton-nucleon interaction potential with a repulsive 
core. We note that the initial scattering wave function 
^ \ is fully antisymmetrized, while in the final chan- 
nel only the bound state wave function ip B is antisym- 
metrized. The transition operator is symmetric relative 
to the nucleon interchange what can be easily seen by 
rewritting it in terms of the kinetic energy operators (see 
Eq. © below). 

Eq. |T]) represents a conventional expression for the 
(d, p) reaction in terms of the volume integral. Recently 
we have developed a new surface-integral based formula- 
tion of the reaction theory with charged particles (see Q 
and references therein) . Transition to the surface integral 
in Eq. ^ can be achieved by replacing the potentials in 
the transition operators by the kinetic energy operators. 
Below we demonstrate it in a few simple and transparent 
steps. First we rewrite Eq. ([1} as 



M dp = <x i f ) MV V B+V B - UpB - V B \^ +} > 



= < 



X }-V B |fr-2*|*i + >>, 



(3) 



(T pB + T B + V pB + V B )y 
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(2) 



where T (T) the total kinetic energy operator acting to 
the right (left) and T= Y,T t . 

i 

Using Eq. ^ we demonstrate that the exact reac- 
tion amplitude Md p is invariant under finite-range uni- 
tary transformation of the wave functions. Such trans- 
formation has been developed in [J] to take into account 
short-range nucleon- nucleon correlations. The problem 



is that in the IPSM, which is used to calculate SFs, the 
nucleus wave function is given by the linear combina- 
tion of the Slater determinants. However, finite num- 
ber of the Slater determinant components cannot repro- 
duce depletion of the nucleus wave function generated by 
a strong short-range repulsive core in the NN interac- 
tion, i.e. the exact wave function contains components 
in the subspace Q beyond of the model space covered 
by the IPSM components. As mentioned above, the uni- 
tary correlation operator method was developed [J] in or- 
der to take into account effectively the repulsive core in 
the iViV-interaction. The correlated wave function in the 
UCOM can be obtained from an uncorrleated state (lin- 
ear combination of Slater determinants) using the unitary 
correlator ^ = C^i>. Here, \P is an uncorrelated state, G 
is the unitary correlation operator (UCO) given by [4[ 



C 



-iG 



where G consists of the radial and tensor 



terms and can be written as G — ^ g{ij)- In this pa- 

i< 3 
per we restrict ourselves only by the two-body operators 
g(ij) assuming that the higher order terms give negligible 
contribution, although the final result is valid for higher 
order correlations as well. Note that the two-body corre- 
lator can generate three- and higher order potential terms 
even if the initial potentials are two-body ones J4|. By 
definition each g(i j) depends on the distance between nu- 
cleons i and j. It is designed to shift away nucleons from 
the classically forbidden region (m 1 — 1.5 fm) generated 
by the repulsive core and disappears in the clasically al- 
lowed region, i.e. g{i j) is a short-range correlator quickly 
going to zero as distance between the correlated nucle- 
ons increases. Similarly three-body correlator disappears 
when one of the three correlated nucleons is moved away 
from the correlated volume. 

Evidently, the application of the UCO conserves the 
norm < \I/|\I/ > = < \IJ| ^ > and the asymptotic normal- 
ization coefficient (ANC), which determines the ampli- 
tude of the tail of the radial overlap function of nuclei A 
and B = (An) (for simplicity we consider the neutron 
ANC): 



lf(r)=<<p A \'PB> « C B An 



r n A 



(4) 



where integration in the matrix element < paIpb > is 
carried over all the internal coordinates of daughter nu- 
cleus A; k is the neutron bound state wave number and 
r nA is the radius connecting n and the center-of-mass of 
nucleus A. Assuming now that both wave functions <pb 
and <pa are correlated we get from Eq. ([4} that 



I A {r nA ) = < Va\^B > = < Va\C a X Cb\&B > 

. iA I - r„ A ->oo B e 

= < <fA\G n A\fB > ~ G At 



TnA 



(5) 



where we used the cluster property of the UCO 

C B =G A G nA . (6) 



Here, Cb is the UCO between nucleons of nucleus B and 
C nA is the UCO between neutron and nucleons of nucleus 
A. 

From the definition of the SF 
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(7) 



we can conclude that the SF, in contrast to the ANC, is 
not invariant under UCO. 

Now using the surface-integral formulation for the re- 
action amplitude we will prove that the reaction ampli- 
tude is invariant under short-range unitary transforma- 
tion of the wave functions. Let \£^ to be the correlated 
wave function in the initial channel (i.e. it contains the 
incident d+A wave), which is related to the uncorrelated 
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Correspondingly the correlated 

channel wave function in the final state $ f = (fB Xf 
can be written in terms of the uncorrelated one as 
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Cb&bXj ) - 



7 - ws 
Now we can rewrite 



(8) 



(9) 
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Note that the distorted wave %V cannot be considered 
as a correlated wave function because it is a solution 
of the Schrodinger equation with the optical potential, 
which depends on the distance between the proton and 
the center-of-mass of nucleus B and doesn't have short- 
range repulsive core. We can distinguish between the out- 
going proton and nucleons belonging to B because there 
is no antisymmetrization between the outgoing proton 
and nucleons of nucleus B, i.e. one can tag the outgo- 
ing proton. Once the transition operator is expressed in 
terms of the difference of the kinetic energy operators we 
can apply the Green's theorem and transform the volume 
integral into the surface one [5| . We rewrite the total ki- 
netic energy operator T as T = T p b+Tb- Since in the bra 
state we have the bound state wave function of nucleus 
B the transformation of the volume integral over internal 
coordinates of nucleus B and extending them to infinity 

leads to the disappearance of the term < | ( T b — T b ) \ >■ 
Another equivalent way to prove that < \(T b — T b)\ > 
disappears is to use the fact that, due of the presence of 
the bound state wave function of nucleus B in the bra 
state, the operator Tb is hermitian. Therefore, taking 
twice the integration by parts we can transform Tb to 



1 b resulting into expression 



M dp =< $ 



[rhfr^-f 



p b-T pB )C p b\W' > 



(+) 



(10) 



Now we split the volume integral (jlQI) into two parts: 
M dp = < $^ } |(^> - f pB )C pB \M +) > 



=Mdp\r pB <R. + Mdp\r pB >R.- 



(11) 



Here t p b is the radius- vector connecting p and the center- 
of-mass of nucleus B. The first term in Eq. (fTTj) describes 



the volume integral, in which r p s < R , while the second 
term determines the contribution of the volume integral 
in which r p B > R. We choose the radius R large enough, 
so that at r p B > R the UCO C p b = 1- The internal 
volume integral we can transform into the surface one 
over the coordinate y p b using the Green's theorem 



L = < /(r)|^ - f\g(r) > \ r < R = -± i? 2 



x / df 



(r)|-r(r)-r(r)|- ff (r) 



\r=R, 



(12) 



where \i is the reduced mass of the interacting particles. 
Taking into account Eq. (]12|) we can rewrite Md p as 



M, 



dp 



M, 



dp 
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,=R 



M, 



dp\r pB >R- 



(13) 



Here, Md p '\ r B =R IS the reaction amplitude in which 
the volume integral over v p b is replaced by the surface 
integral taken over $7 r B at r p s = R. The integrations 
over all other independent coordinates are carried over 
without limitations. The operator C p b is expressed in 
terms of the two-body operators g(pj), which connect 
the proton and nucleon j G B. Each g(pj) consists of the 
radial and tensor parts, which depend on r pj - = t p b — i'j- 
Here, Tj i s the radius- vector connecting nucleon j and 
the center-of-mass of nucleus B. Since in the final state 
nucleus B is in a bound state, the integration over all 
Tj, j G B, is limited. Hence at large enough t p b radius 
r p j will also be large enough to exceed the correlation 
radius, i.e. g{pj) = 1. Thus at large enough R the UCO 
C p b can be replaced by unity, i.e. Md p ' \ rpB =R becomes 
insensitive to the UCOs. The second term, Md p \ rpB >R, 
is the reaction amplitude in which the volume integral 
over y p b is taken from r p s = R and all the integrations 
over other independent variables are carried out over the 
whole configurational space. Thus starting from the tra- 
ditional expression (TTJ) and using the surface integration 
we have proved that the exact {d,p) reaction amplitude 
is invariant under the finite-range unitary correlations. 

The invariance of the reaction amplitude is true only 
if the exact formulation is used, i.e. the exact wave func- 
tion in the initial state (the post-form) or the exact final 
state wave function (in the prior- form). The asymptotic 
behavior of the exact wave function is given by the inci- 
dent wave plus the elastic, inelastic, rearrangement and 
breakup scattered waves [6| . The amplitude of each scat- 
tered wave is the reaction amplitude for transition initial 
channel i —¥ the final channel corresponding to this scat- 
tered wave. For the rearrangement channel A{d,p)B the 
amplitude of the outgoing wave in the exit channel is the 
reaction amplitude given by Eq. (TTJ). The knowledge of 
the exact wave function for the collision d + A assumes 
that we know its asymptotic behavior in all the asymp- 
totic regions Q, i.e. we know the amplitudes of the all 
opened rearrangement and breakup channels. 

The invariance of the reaction amplitudes under the 
UCOs is understandable: these amplitudes determine 



the asymptotic behavior of the scattering wave function 
which is not affected by the finite range UCOs. Although 
we have proved the invariance of the asymptotic behav- 
ior of the wave function in the two-fragment channels 
below we prove that it is also the case for the breakup 
reactions. Thus, the exact reaction amplitude cannot be 
used to determine unambiguously the SF which is not 
invariant under UCO. Moreover, the exact reaction am- 
plitude is not parametrized in terms of the SF, which 
can be calculated as the square of the norm of the over- 
lap function or using the method discussed in Q- In 
contrast to the exact approach, the conventional DWBA 
is not invariant under UCO. The DWBA is designed to 
determine the phenomenological SFs by comparing the 
calculated differential cross section with the experimen- 
tal one. To express the DWBA amplitude in terms of the 
SF drastic approximations are done for the initial exact 
scattering wave function and the transition operator V p b ■ 
In particular, in the initial scattering wave function only 
the incident wave is left, while all other opened channels 
are neglected and coupling to them is taken into account 
only effectively in terms of the initial distorted wave. The 
approximation of the initial scattering wave function by 

the channel wave function <pd ipA x\ hi Eq. (|T0|) leads 
to zero for the reaction amplitude, in which transition 
operators are written in terms of the kinetic energy op- 
erators: once we write the transition operators in terms 
of the kinetic energy operators, the volume integral can 
be transformed in to the surface one and the radius of the 
surface should be large enough to encircle the volume in 
which reaction occurs. To calculate the surface integral 
the correct asymptotic behavior of the exact scattering 
wave function is required. Moreover as the radius of the 
surface goes to infinity (for the (d,p) reaction it would be 
r p B = R — > oo ) only the contribution from the outgoing 
wave in the exit channel p + B will survive for the (d, p) 
reaction but this wave is missing in the DWBA. 

Let us write now the DWBA reaction amplitude (in 
the post-form) taking into account UCOs: 



M. 



DW 



< $ { f> C B l \V pB - U pB )\C pn C A *W >, (14) 



>}->/" 



where we took into account that the channel wave 



function in the initial state $ 



(+) 



,(+) 



C pn C A & + \ 



i(+) - 



,(+) 



= fd <PA Xi 

^f' = (pd&AXV • I n tue adiabatic 
DWBA (ADWBA) V pB - U pB is replaced by V pn . Then 
we get the ADWBA amplitude taking into account 
UCOs: 



Mg w =< iW-j V pn C pn \^ > 



=< 






X) ' lA\ V pnC pn \(fdX. 



(+) 



> 



(15) 



Since the (d, p) reactions are dominantly peripheral 
(about 70 — 80% contribution to the reaction amplitude 
comes from the nuclear exterior), the probability to find 
the transferred neutron close to nucleons in A at dis- 
tances of the range of the UCO is negligible Q, i.e. the 
effect of the correlator C~\ on the DWBA amplitude can 



be neglected. The main impact of the UCO is replace- 
ment of V pn by V pn C pn . However, this effect is too small 
to significantly affect the DWBA amplitude. 

Breakup (d,pn) reaction amplitude. Here we demon- 
strate that the deuteron breakup reaction amplitude is 
also invariant under UCO transformation. Usually to 
consider the breakup amplitude one starts from the prior- 
form: 



,(-),.(-) 



,-,(-) „(-) ,~„ 



v(-) 



M d pn = < */ \VdA - UdA\<Pd¥>AX 



(+) 



>=< ¥ f ] \lf 






T~ + 7 ?\<S> ( ; +) > 



(16) 



Here, evidently, ^ \ is the exact scattering wave func- 
tion in the final state, which has the three-body incident 



wave p + n + A and H *$> 
tonian of the system p 



(-) 

/ 

n 



EV 
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/ 



H is the Hamil- 



A. The total kinetic energy 



operator can be written as T — T pn + TdA + Ta, where 
T pn and Ta are hermitian because in the initial state 
deuteron and A are bound. Hence, 

M dpn = < ¥f ] \ - T~ dA + f dA \(pd(PAX4 +) > 

= < ^C- 1 ] - f dA + TdAlCdCAtfdVAX^ > 

= < */ _) c7il - ^d A + ?dA\<Pd<PAjci +) > ■ (17) 

where C = C A C d C d A, C d = C pn . The integral over T d A 
is not hermitian because the relative d — A motion in 
the initial and final state is unbound. Note that we can 
tag the proton and neutron belonging to the deuteron in 
the initial state because the wave function in the initial 
state is antisymmetrized only over exchange of nucleons 
belonging to nucleus A. 

We can repeat now the steps used when considering the 
(d,p) amplitude. We divide the volume integral over TdA 
into the internal and external parts, and transform the in- 
ternal volume integral into the surface one. If the radius 
of this surface is large enough the UCO CdA — C p a C n A 
can be replaced by the unit operator. Evidently the ex- 
ternal volume integral in such a case also doesn't depend 
on CdA- Thus we proved that the deuteron breakup am- 
plitude is also invariant under the UCO transformation. 
Note that we can include the three-body correlator C pn A 
into CdA and it still turns into unity when TdA becomes 
large enough. Since we have proved before the invari- 
ance of the reaction amplitudes for the rearrangement 
channels we can conclude that asymptotic behavior of 
the exact $> ' in general is invariant under the UCO 
transformation. We note that for the nucleon knockout 
in heavy-ion collisions the invariance of the reaction am- 
plitude under the UCOs can be proved in a similar way. 

The DWBA amplitude can be obtained from Eq. (fTo) 

by replacing the exact wave function \Jr. by the three- 
body channel wave function, which is , in contrast to the 
two-body case, is not determined uniquely and is dif- 
ferent in the different asymptotic regions. If we take 



the channel wave function in the final state as $ 



(-) 
/ 



XW XdA VA = C pn C A X P n X d A <f 'A (or its CDCC ex- 
tension) and take into account that VdA in the DWBA is 
replaced by the interaction potential, which depends on 
the distance between the center-of-mass of d and A, then 
we see immediately that the DWBA amplitude doesn't 
depend on the UCOs. Thus the DWBA amplitude for the 
(d,pn) reaction is invariant under UCO transformation. 
(e, e' p) reaction amplitude. — The B(e, e' p)A process 
amplitude in the prior-form is given by (similar to the 
(d,pn) reactions) 



M ee , p = <& f ) \V eB - U eB \ VB xi +) > 



= < ¥ f ] \ -1r + '?| 



^x! +) > 



(18) 



Note that this expression is valid even if the three-body 
potentials V ep j are included in addition to the two-body 
potential V ep . Introducing the correlators we get 



M ee , p = < ¥~ } C-\C A l I - T~ + f\C pA C A £bX- +) > 

(19) 



/ ^pA^A 
1 eB 



<¥ f ^-^TeB+^c 



B|^i?xi +) >) 



,(+) 



where T = T e b + TpA + Ta and \\ is the electron 
distorted wave in the initial channel. To obtain the last 
equation we took into account that operators T p a and 
Ta are hermitian because B is in the bound state in the 
initial channel. Thus the exact (e, e'p) amplitude is in- 
variant relative UCOs. The DWBA will be considered 
els where. 

In conclusion, we have demonstrated that the 
(d,p), (d,pn) and (e, e'p) reaction amplitudes determin- 
ing the amplitudes of the asymptotic terms of the exact 
scattering wave functions in the corresponding channels 
are invariant under finite range UCOs while the SFs are 
not. Moreover, the exact reaction amplitudes cannot be 
parametrized in terms of the SFs, which should be cal- 
culated from the overlap functions rather than from the 
reaction amplitudes. Hence, the exact reaction ampli- 
tudes cannot provide a tool to determine the SFs. Only 
drastic simplifications of the reaction amplitudes lead to 
the DWBA amplitudes, which are parametrized in terms 
of the " SFs" determined by comparison with the experi- 
mental data. However, in reality, these "SFs" are nothing 
else but proportionality coefficients between the DWBA 
cross sections and experimental data. The DWBA am- 
plitudes for the (d, p) reactions are sensitive to the finite 
range UCOs but, due to the peripheral character of the 
(d,p) reactions, this sensitivity is suppressed, while the 
DWBA amplitudes for the (d, pn) reactions are invariant 
under the UCOs. 
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